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Abstract 

The communication cost of simulating probability distributions 
obtained by measuring quantum states is a natural way to quantify 
quantum non-locality. While much is known in the case of bipartite 
entanglement, little has been done in the multipartite setting. In this 
paper, we focus on the GHZ state. Specifically, equatorial measure- 
ments lead to correlations similar to the ones obtained with Bell states. 
We give a protocol to simulate these measurements on the n-partite 
GHZ state using 0(n 2 ) bits of communication on average. 

1 Introduction 

The issue of non-locality in quantum physics was raised in 1935 by Einstein, 
Podolsky and Rosen [5]. Thirty years later, John Bell proved that quantum 
physics yields correlations that cannot be reproduced by classical local hid- 
den variable theories [2]. This momentous discovery led to the more general 
question of quantifying quantum non-locality. Not only is this question rele- 
vant for the foundations of quantum physics, but it is directly related to our 
understanding of the computational power of quantum resources. 

A natural quantitative approach to non-locality is to study the amount of 
resources required to reproduce probabilities obtained by measuring quantum 
states. In this paper, we consider the simulation of these distributions using 
classical communication. This approach was introduced independently by 
several authors [HI EJ HO]. It led to a series of results, culminating with 
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the protocol of Toner and Bacon to simulate von Neumann measurements 
on Bell states with a single bit of communication [TT]. Later, Regev and 
Toner extended this result by giving a simulation of binary von Neumann 
measurements on arbitrary bipartite states using two classical bits [9]. 

We focus here on multipartite entanglement, and more specifically on 
GHZ states [7j. Unlike the bipartite case, which has been the topic of inten- 
sive investigation, the simulation of multipartite entanglement is still teeming 
with major open problems. 

The easiest situation arises in the case of equatorial measurements on 
a GHZ state because all the marginal probability distributions obtained by 
tracing out one or more of the parties are uniform. Hence, it suffices in this 
case to simulate the n-partite correlation, henceforth called the full corre- 
lation. (Once this has been achieved, all the marginals can easily be made 
uniform [6].) Making the best of this observation, Bancal, Branciard and 
Gisin have given a protocol to simulate equatorial measurements on the tri- 
partite and fourpartite GHZ states at an expected cost of 10 and 20 bits 
of communication [TJ. However, the amount of communication entailed by 
their protocol is unbounded in the worst case. More recently, Branciard 
and Gisin impoved this in the tripartite case with a protocol using 3 bits of 
communication in the worst case [I]. 

In this paper, we give a protocol to simulate equatorial measurements 
on the n-partite GHZ state. For any measurements, our protocol has an 
expected cost of 0(n 2 ) bits of communication, where the expectation is taken 
over the inner randomness of the protocol. 

The paper is organized as follows. In the Section we give the structure 
of the distribution arising from equatorial measurements on GHZ states, and 
then show how to simulate it. The main technical tool that we use is a 
protocol to sample uniform vectors on connected subsets of the circle. We 
call this task Uniform Vector Sampling. The protocol to sample those vectors 
is given in Section |3j 

The structure of our simulation is inspired by the protocol proposed by 
Toner and Bacon to simulate von Neumann measurements on Bell states. 
It can be divided in two parts. In the first part, the players communicate 
to sample shared random vectors on the sphere. The distribution of these 
vectors depend on the player's input. In the second part, they apply a post- 
processing to compute their outputs. While the sampling part in the protocol 
of Toner and Bacon is simple, it is more involved in our case. The goal of 
Uniform Vector Sampling is to sample vectors with the correct distribution. 
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2 Simulating equatorial measurements 



We consider the family of GHZ states |*„) = ^(|0 n ) + |l n )), and the dis- 
tribution generated by the following process: n players each receive a qubit 
of l^n). Each player applies a bipartite measurement to its share. Let 
Oj G {—1,1} denote the output of the i th player. The problem is to simulate 
the probability distribution over the player's output using hidden variables 
and communication. 

The measurement operators corresponding to equatorial measurements 
are on the equator of the Bloch sphere and therefore, can be parametrized 
by a single polar angle. Denote oij the angle corresponding to the i th player's 
measurement. It is known that the distribution arising from such measure- 
ments is fully characterized by the full correlation (see e.g. pQ). 

Proposition 1. The distributions of the outputs {o{\ is characterized by the 
following relations: 

• The full correlation is given by E \\Xi=i °i] = cos (Sr=i a i)- 

• The marginal distributions are given by E [ILes 0i ] = ® f or a ^ $ £ M • 

We denote S 2 the sphere in dimension 3 and S 1 the sphere in dimension 
2. A vector on S 1 is parametrized by a single polar angle, or equivalently a 
real number modulo 2tt. An interval on S 1 is a connected subset of S 1 or 
equivalently of M/27rZ. 

Our simulation is based on a procedure to sample uniform vectors on 
intervals of S 1 , when the description of this subset is shared among sev- 
eral players. For k,n > 1, we introduce the following task, called Uniform 
Vector Sampling and denoted UVS(n, k). The n players each receive the 
angles «i, . . . , a n , respectively. Each player computes a message depending 
on his input and on a public random variable r and sends it to a referee. 
At the end, the referee has to output a uniform angle 6 on the interval 
EiLi a i ~ 71 /2 k , Yli=i a i + 7r /2 fc ] • We measure the communication cost of a 
protocol for UVS by considering the total length of all messages sent from 
the players to the referee. 

Theorem 1. For any n,k > 1, there exists a protocol for UVS(n, k) with 
expected communication cost at most n(n + k). 

We now show how to simulate equatorial measurement on GHZ states, 
given a protocol for UVS. Toner and Bacon proposed a simulation of binary 
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measurements on Bell states, using a single bit of communication. In the 
bipartite, the correlation between the player's output is a scalar product of 
two vectors on S 2 . We sketch their protocol. For a vector a 6 S 2 , denote 
S + (a) the half sphere centered on a, S + (a) = {A G S 2 : (a, A)}. We denote 
the sign function with range { — 1, +1} by sgn. Toner and Bacon prove the 
following theorem. 

Theorem 2 ([11]). Let a,b be vectors in S 2 , and Ai,A 2 be uniformly dis- 
tributed on S + (a). Then 

E[sgn(A 1 + A 2 ,6)] = (a, b). 

To complete the simulation, it suffices to notice that shared uniform vec- 
tors on S + (a) can be sampled efficiently by players using shared randomness 
and communication, even if only one player has a full description of a. The 
idea is to first sample a uniform random vector on the sphere, and then the 
player that knows a tells the other if he has to flip the random vector in order 
to get a vector in S + (a). This requires to send a single bit of communication. 

Our simulation is based on the following observation. Consider d = 
(cos Y17=i a i> s * n J2^=i a i> 0) an d a n = (cos a n , — sin a n , 0). These are unit 
vectors on S 1 , embeded in M. 3 to apply Theorem |2j For these vectors, we 
have (d,a n ) = cos^™ =1 a;. Therefore, if Ai and A2 are two vectors sampled 
uniformly on S + (d), Theorem [2] gives E[sgn(Ai + A 2 , a„)] = cos^c^. 

We now describe the simulation in more details. The players are denoted 
Ax, . . . , A n . Before receiving their inputs, they prepare a shared variable r, 
used for UVS. In addition, they prepare some shared uniform random bits 
bi G { — 1, +1} for i — 1, . . . , n — 1. In our simulation, we only need to apply 
Uniform Vector Sampling with k = 1. 

1. For i — 1, . . . , n — 1, the players Ai run the protocol for UVS(n — 1, 1), 
sending their messages to A n . 

2. Using the messages he received, A n sets 9\ uniform on the interval 

Er=i 1 ^-^/2,Er=i 1 ^+^/2]. 

3. The players repeat steps [1] and |2] to allow A n to sample another angle 
#2 with the same distribution. 

4. Player A n samples U\ and U2 uniformly on [—1, 1] and for i = 1,2, sets 

ifi = arccos Ui 

Xi = (cos#; cos^j, sin6*i cos (fx, sin (pi). 
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5. For % = 1, . . . , n — 1, the player A { outputs Oj = 6j. 

6. The player A n outputs o n = (Yl™=i h) ■ sgn(Ai + \2,a n ), where we 
defined a n = (cosa„, — sina n , 0). 

After step EJ player A n has the complete description of two angles 9i 
and #2 uniformly distributed on ai ~ f > ^i=i ai + f ] • The P ur P ose 

of step H] is to transform the angles in uniform random vectors on S + (d), 
where d is the vector with coordinates (cos Y^=i a *> sm Yll=i oii,0). Since d 
is on the equator, it is sufficient to assign a random latitude to the vectors 
whose longitudes are 6\ and 9 2 . Finally, after steps [3] and El we have 

E Y\ Oi = for any S £1 [n] , 

n 

EjJ 0i = E[sgn(Ai + A 2 ,a„)] 
i=i 

n 

= cos ctj by Theorem [2l 

i=i 

Sampling the angles d\ and # 2 , can be done with 0{n 2 ) expected bits of com- 
munication. Therefore, the whole protocol can be done with 0(n 2 ) expected 
bits of communication. 

3 Uniform Vector Sampling 

The goal of this section is to prove Theorem[TJ Observe that in the simulation, 
we only need the case k — 1. Nevertheless, our inductive proof requires to 
prove the stronger statement given in Section [21 

The base case: n = 1 

For n = 1, there is a single input a,\. Fix k > 1. Let S be chose uniformly at 
random on S 1 . The player sends 

t = mm{ieN:6 + i^ =I e + 

to the referee, who computes 6 = 5 + The resulting angle 9 is uniformly 

distributed on [a\ — ot\ + Jj]. Notice that since t G [2 fe — 1], the length of 
the message is at most k. 
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The induction step 

Let n > 1. The following Lemma is the main technical tool that we use for 
the induction. It explains how to generate uniformly distributed variables 
from specific non-uniform ones. We first give and prove Lemma [U and then 
use it to prove the induction. 

Lemma 1. LetT>^ denote the uniform distribution on [0,1/2*] andT>^ de- 
note the uniform distribution on [1 — 1/2 4 , 1]. Let T> be the distribution on t 
defined by the following procedure: 

• Pick an integer i > with probability l/2 l+1 , andr uniform in{— 1, +1}. 

• If r = —1, sample ti,t 2 ~ T>^ . 

• Otherwise, sample ti,t 2 ~ T>f . 

• Set t = t\ + t 2 . 

Then T> is the uniform distribution on [0, 1]. 

Proof Denote Ur = [0, 1/2*] and Uf = [1 - 1/2% 1]. We define the density 
functions associated to the distributions T>f and T>^ , 



where * denotes the convolution product of two functions. By direct calcu- 
lation, we have 




By definition, the density pi of t\ + t 2 for a fixed i is 



Pi = \tfi*fi +fi *fi )> 




22(1+1)3. if x e [o,i/2 i+1 ], 

2 l+2 - 2 2 ^x if x G [l/2 i+1 , 1/2*], 



otherwise. 



and 



(ft*ft){x) = (fr*fr)0--*)- 
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Let p denote the density of the distribution V. We now calculate p(x). Notice 
that fo = fo, and for % > 0, /~ and have disjoint supports. Assume that 
x < 1/2 (the other case is similar). In that case, f^{x) = for any % > 0. 
Let j = max{/ : x G [0, 1/2 J ' ]} and notice that ff(x) = for any i > j. We 
have 



The induction hypothesis is that for any k > 1, it is possible for n — 1 
players to each send a message to another party such that he outputs an 
angle 9 uniformly distributed on [X^] 1 a « ~~ K /2 k , YH=i a i + 7r /2 fe ]- 

Before receiving their inputs, the players prepare the following random 
elements: 

• an integer j > chosen with probability p(j) = 1/2 J+1 , 

• b uniform in { — 1 , + 1 } , 

• the random elements required to run UVS(n — 1, k + j + 1), 

• the random elements required to run UVS(1, k + j + 1). 

The protocol proceeds as follows. The n— 1 players first players send to the 
referee the messages corresponding to UVS(n— 1, k+j + 1). The referee uses 



them to prepare $i uniform on EILi* cti — Tx/2 k+j+l , Y17=i a i+ 7r /2 k+j+1 ]. The 




\ ■ \(fo * /o + )(^) + 2^r(/< * fi ) 

i=0 

x + - Vi2 2 (' +1 )x + l(2^-2Wi) 
L «=o 




x + (2 j - l)x + 1 - 2 j x 
1 



which concludes the proof. 



□ 



7 



0.6 



0.4 



0.8 



0.2 




0.2 



0.4 



0.6 



0.8 



1 - 



0.2 0.4 0.6 0.8 



Figure 1: The density functions and /~, for i < 3 and i < 8. Each density 
function is scaled down by the probability of sampling it in Lemma [TJ /q" 
and Jq are equal. The dashed curves represent the sum of the represented 
density functions. 

n-th player sends to the referee the message corresponding to UVS(1, k + j + 
1). The referee uses it to prepare 9 2 uniform on [a n — 7r/2 fc+J+1 , a n +7r/2 fe+J+1 ]. 
Finally, the referee outputs 9 = 9\ + 9 2 + b ■ ^ (l — . 

We now analyze the protocol and prove that 9 is uniform on the interval 
[X)iLi a i ~~ 7r /2 fc , Y^=i a i + 7r/2 fe ] . To apply Lemma [H we need to rescale 
the random variables 9\ and 9 2 . We split the term b ■ ^(1 — in two parts 
and think of each as a shift of 9\ and 9 2 in a direction that depends on the 
bit b. Each angle is shifted in the same direction. 

Let v^j and v^j be uniform random variables on [0, 1/2- 7 '] and and Uj- 
be uniform random variables on [1 — 1/2 3 , 1]. Let Tij denote the random 
variable 9\\j, that is, the random variable generated by UVS(n — l,k + j + l) 
for a fixed value of j. The shifted random variable Tij + fe^rr (l ~~ ^j) is 



• either on [Y^=i «i ~ 7r/2 fe+1 , Yn=i a i ~ ^/2 k+1 + n/2 k+ i] if b = -1, 

• or on [EL"/ «i + ^/2 fc+1 - 7r/2* +i , E"=i <*i + ^/2 fc+1 ] if 6 = +1. 



uniform 



Using, uj". and t> x we can rewrite 




E^ 1 «i - vr/2 fe+1 + ^ • rc/2 k if 6 = -1 
££"1 «i - vr/2 fe+1 + u+ ■ n/2 k if 6 = +1 
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Similarly, let T 2j - denote the random variable 02\j- Using v 2 „• and v^j, it can 
be written 

T I h I (\ M - { a n-n/2 k + l + v^, r n/2 k if 6 = -1 
2J 2 fc +! V 2>J ]a n - 7r/2 fc+1 + «+ • 7t/2 k if 6 = +1 

For the sum, we get the expression 

9 = Ty + T 2J + 6 • £ (l - 1) = £ a, - n/2 k + v hb • Tr/2^ 1 

^ ' i=l 

where is the sum of vfj and v£j if b = +1 and the sum of v^j and 
v^j if 6 = —1. According to Lemma [U when taking the expectation over j 
and b, Vjt> is uniform on [0, 1]. In consequence, 9 is uniform on the interval 

It remains to bound the expected length of messages. Denote l n ^ the 
expected sum of the messages length. We already know that < k for 
any k. Fix n > 1. Analyzing our protocol, we get the induction: 

ln,k = 7^TT^n-l,k+j+l + k,k+j+l), 



2? 

j>o 

i>o j>o 
j>o 



The induction hypothesis is that l n -i,k+j+\ < (n — l)(n + k + j). We plus 
this expression and get 

i n +i,k < 2^1 #n J+ fc ' 

< (n - l)(n + k + l) + k, 

< n(n + fc), 

which concludes the proof. 
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4 Conclusion 



We gave a protocol to simulate equatorial measurements on the n-partite 
GHZ state, using 0(n 2 ) bits on average. Our protocol is in two parts. Firstly, 
we reduce the problem to sampling vectors on regions of the S 1 . Secondly, 
we give a procedure to sample the vectors, called Uniform Vector Sampling. 
This scheme is inspired by the protocol of Toner and Bacon to simulate von 
Neumann measurements on Bell States. 

Our work leads to an obvious question. Is it possible to transform our 
protocol into a protocol that is bounded in the worst case? To solve this 
question, it enough to give a protocol for UVS that use bounded communi- 
cation in the worst case. Uniform Vector Sampling could also be considered 
as a task of independent interest or be applied in other contexts. 

Our work, like others on the same topic, considers only equatorial mea- 
surements. The simulation of more general measurements is an intrigu- 
ingly hard question. The main difference is that they lead to non-uniform 
marginals. In the bipartite case, an analogous problem arises when consid- 
ering non-maximally entangled states. It may seem that modifying local 
marginals is easy once the correlation is simulated. Unfortunately, local 
transforms usually also modify the full correlation. 
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